VERY SINGULAR SOLUTIONS 
FOR THIN FILM EQUATIONS WITH ABSORPTION 



V.A. GALAKTIONOV 

Abstract. The large-time behaviour of weak nonnegative and sign changing solutions 
of the thin film equation (TFE) with absorption 

m = -v- (M n vAu) - M P "V 

where n 6 (0, 3) and the absorption exponent p belongs to the subcritical range 
pe (ra + l,i> ), with p = 1 + n + ±, 

is studied. Firstly, the standard free-boundary problem with zero-height, zero contact 
angle and zero-flux conditions at the interface and bounded compactly supported initial 
data is considered. Very singular similarity solutions (VSSs) have the form 

u s (x,t)=r^f(y), y = x/iP, = > 0. 

Here / solves the quasilinear degenerate elliptic equation 

-v • d/rvA/) + /3v/ • y + ^ f - i/r 1 / = o 

that becomes an ODE for N = 1 or in the radial setting in M. . By a combination of 
analytical, asymptotic, and numerical methods, existence of various branches of similarity 
profiles / parameterized by p is established. Secondly, in parallel, changing sign VSSs of 
the Cauchy problem are described. 

This study is motivated by the detailed VSS results for the second-order porous 
medium equation with critical absorption (u > 0) 

u t = V • (u"Vu) - uP ia K" x K + , p=l + n+^,n>0, 

which have been known since the 1980s. 



1. Introduction: very singular similarity solutions 

1.1. On the PME with critical absorption. Second-order quasilinear diffusion equa- 
tions with absorption are typical for combustion theory. The most well-known model, 
which became a canonical object of intensive investigation in the 1970-90s, is the porous 
medium equation (PME) with absorption 

(1.1) u t = V • (u n Vu) - u p in R N x R + (u > 0), 

where n > and p G K are given parameters. Mathematical theory of such PDEs 
was founded by A.S. Kalashnikov in the 1970; see his fundamental survey [30] for the 
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full PME history in the 1950-80s. Besides new phenomena of localization and interface 
propagation, for more than twenty years, the PME with absorption (II. II) was a crucial 
model for determining various asymptotic patterns, which can occur for large times t ^> 1 
or close to finite-time extinction as t — > T~ (for p < 1). For ( 11 . lit . there are several 
parameter ranges with different asymptotics of solutions such as 

p>Po = l + n + jj, p = Po, l + n<p<p , p=l + n, 
1 < p < 1 + n, p = 1, 1 — n < p < 1, p = 1 — n, p < 1 — n, 

etc.; see details, references, and history in the book [221 Ch. 5, 6]. A transitional behaviour 
for (11.11) occurs precisely at the first critical absorption exponent 

(1.2) p = l + n + %, 

where logarithmically perturbed asymptotic patterns occur. VSSs, as special key solu- 
tions, are known to appear in the subcritical range p e (n + l,po), while for p > p Q , 
the asymptotics of solutions as t — > oo correspond to pure diffusive PME, where the 
absorption term becomes negligible. 

We will consider a new quasilinear parabolic model obtained by adding to the standard 
thin film operator an extra absorption term. This creates a non-conservative evolution 
PDE with the differential operator that is only "partially" divergent. We are going to 
show that, nevertheless, this higher-order nonlinear PDE exhibits several evolution fea- 
tures that are surprisingly similar to those for the PME with absorption (11.11) . Moreover, 
we also claim that mathematical features and properties of both equations are also similar 
but, indeed, the results are much more difficult to justify for the higher-order case. Sev- 
eral conclusions of the present paper remain interesting and difficult open mathematical 
problems, which we would like to stress upon. 

We claim that many of the presented results and ideas cannot and will not be proved 
rigorously in a few years to come, but a reasonably detailed qualitative understanding of 
such popular nowadays models and nonlinear PDEs is essential. 

1.2. The TFE with absorption. Out goal is to show that VSS phenomena exist in the 
new model that is the thin film equation (TFE) with absorption 

(1.3) u t = -V • (\u\ n VAu) - \u\ p -\ 

where n > and p > 1 are fixed exponents. For convenience, it is written for solutions 
of changing sign to be also studied. Extra absorption terms in such TFEs are to model 
the effects of evaporation or permeability of the surface. We refer to [35] and [Tj, [27J for 
derivation of various TFE models including the non- conservative cases. Other related 
details and references on such PDEs can be found in the previous paper [21 J on the TFE 
fll.3p devoted to the delicate case of the critical absorption exponent 

(1.4) p Q = l + n+±. 

In the present paper we study equation (11.31) in the complementary subcritical range 

(1.5) n + l<p<p = l + n + j ! . 

2 



All necessary key references on various results on modern TFE theory that are necessary 
for justifying principal regularity and other assumptions on solutions will be presented 
below and are partially available in [21] . See also [231 Ch. 3] , where further references are 
given and several properties of TFEs with absorption are discussed. 

We will consider two main problems for the TFE ( 11. 3ft : 

(I) the standard free-boundary problem (FBP) admitting nonnegative solutions, and 

(II) the Cauchy problem with another functional setting and solutions of changing sign. 

(I) For the FBP, equation (I1.3P is equipped with zero-height, zero contact angle, and 
zero-flux (conservation of mass) free-boundary conditions 

(1.6) u = Vit = v • (m"VAm) = 

at the singularity surface (interface) T [u], which is the lateral boundary of suppii with 
the outward unit normal v. Bounded, smooth, and compactly supported initial data 

(1.7) u(x, 0) = u (x) in T [u] n {t = 0} 

are added to complete a suitable functional setting of the FBP. We then assume that 
these three free-boundary conditions give a correctly specified problem for the fourth- 
order parabolic equation, at least for sufficiently smooth and bell-shaped initial data, 
e.g., in the radial setting; see references below. 

(II) For the Cauchy problem (CP), we need solutions exhibiting the maximal regular- 
ity at interfaces, and this demands oscillatory character of such solutions. We refer to the 
book [231 Ch. 3] and [T5l [T7] for details concerning definitions and details on oscillatory 
changing sign solutions of various TFEs and other nonlinear higher-order PDEs. 

1.3. On the pure TFE and similarity solutions. In what follows we will need well- 
known similarity solutions of the usual divergent TFE with mass conservation 

(1.8) u t = -V- (\u\ n VAu). 

Earlier references on derivation of such fourth-order TFE and related models can be found 
in [2S1 EH], where first analysis of some self-similar solutions for n = 1 was performed. 
Source-type similarity solutions of (11. 8p for arbitrary n were studied in [6] for N = 1 and 
in [18] for the equation in Mr. More information on similarity and other solutions can be 
found in [UEICF]. TFEs admit non-negative solutions constructed by "singular" parabolic 
approximations of the degenerate nonlinear coefficients. We refer to the pioneering paper 
[2] and to various extensions in [131 HU EH [38] and the references therein. See also [29] for 
mathematics of solutions of the FBP and CP of changing sign (such a class of solutions 
of the CP will be considered later on). 

In both the FBP and the CP, the source-type solutions of the TFE ( II. 8p take the form 

(1.9) u s {x,t) = t~ m F(y), y = x/t 13 , with exponent (3 = 

where, in the FBP, F(y) > is a radially symmetric compactly supported solution of the 
ODE PUB] 

(1.10) B(F) = -V • (F n VAF) + [3VF ■ y + (3NF = 0. 
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In the case n — 1, the similarity profile for the FBP is given explicitly, 



(1-11) F{y) = c {a 2 -\y\ 2 ) 2 , c = 8(jV+2 ) (JV+4) , a > 0, 

which was first constructed in [36] . Most advanced results of asymptotic stability for the 
TFE (11.81) were obtained for the similarity solutions with profiles (II. lip for n = 1, where 
the rescaled equation is a gradient system. See a full account of such studies in |ll] (L 1 - 
convergence) , [10] (if ^convergence) and also |25j for more general properties including 
existence and uniqueness for the FBP (cf. [T5J § 6.2], where these questions were treated 
in the von Mises variables). For other values of n ^ 1, the results are much weaker 
and are not complete. This emphasizes how difficult the non fully divergent TFEs with 
non-monotone operators are for the qualitative asymptotic study. 

For the Cauchy problem, there exists a unique (up mass-scaling) oscillatory similarity 
profile of (ll.lOp for not that large n G (0, n^), = 1.759... [15], while for n G (0,1), 
existence and uniqueness are straightforward consequences of the result in [5] on oscilla- 
tory solutions. See also some details in [231 § 3.7] and interesting oscillatory features of 
similarity dipole solutions of the TFE (11.81) observed in [8]. Stability of such sign chang- 
ing similarity solutions is quite plausible but was not proved rigorously being an open 
challenging problem. 

1.4. On main results: comparison of logarithmically perturbed asymptotics 
for p = po and VSSs for p < p . The critical case p = Po was studied in [2T], where 
it was shown that the asymptotic behaviour of solutions is governed by a In t-perturbed 
similarity solution (11.91) . For instance, for n = 1, relying on the explicit representation 
(11.111) and good spectral properties of the corresponding self-adjoint linearised rescaled 
operator, it was shown that, for p = p Q = 2 + the TFE with absorption (II. 3p admits 
asymptotic patterns of the following form as t — > oo: 

(1.12) u(x,t) ~ (tlnty^F^x/t^lnt)-^ 4 ) (/3 = ^) , 

where F* is a fixed rescaled profile from the family (11.111) with a uniquely chosen parameter 
a = a* (AT) > 0. For the semilinear case n = 0, i.e., for the fourth-order parabolic equation 
written for solutions of changing sign 

(1.13) u t = -A 2 u-\u\ p - 1 u, 

the critical behaviour like (I1.12p is known to occur at the critical exponent p = 1 + 4 
[19], which is precisely (11.41) with n = 0. In this case, the centre manifold analysis also 
uses spectral properties of a non self-adjoint linear operator studied in [T2j . 

It turns out that, in the subcritical range ( ll.5p . the generic behaviour of compactly 
supported solutions needs a different similarity description leading to the VSSs. In the 
next Section [21 we describe the VSSs in the subcritical range p G (n + l,Po) f° r the 
FBP. In the final Section [31 we discuss VSS structures for the Cauchy problem admitting 
maximal regularity solutions of changing sign. In general, we show the branching of the 
VSS profiles from the similarity patterns F of the pure TFE. 
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2. VSSS FOR THE FBP IN THE SUBCRITICAL RANGE p G (n + l,p ) 



In the range (11.5ft . the TFE with absorption (11. 3ft admits the self-similar very singular 
solution (VSS) of the standard form 

(2.1) u s (x,t)=r^f(y), y = x/t?, /? = ^±i) > , 
where / solves the nonlinear elliptic equation 

(2.2) -V • (|/rVA/) + 0V/ . „ + -L. / - = 0. 

For n = and p E (1, 1 + 4), the existence of a finite number of oscillatory profiles / 
in M N corresponding to the Cauchy problem for (j 1.13p was detected in [23] by using a 
p-bifurcation analysis. 

For the current FBP, the solvability of (12.21) and existence of a nonnegative compactly 
supported / is unknown even in radial geometry, and we are going to present some analytic 
and numerical evidence concerning this. In ID, (12.21) is an ODE, 

(2.3) -(\f\nf»y + f]fy + -L 1 f- |/r7 = foryG(0,y ), f(0) = / w (0) = 0, 

where we put two symmetry boundary conditions at the origin. At the interface y — yo, 
the solution expansion is as follows (see [6l fl5l [18]): 

(2-4) f(y) = C (yo - y? + p-^A-an) fa> " v)^ 1 + °«)> 

where yo > and Co > are arbitrary two parameters. Therefore, in general, matching 
the 2D bundle (12.41) . comprising two positive parameters yo and Co, with two symmetry 
conditions in (12.31) cannot give more than a countable number of similarity FBP profiles 
{fi} (provided the parameter dependence is analytic). 

Figure [T] shows typical first positive symmetric VSS profiles fo(y) constructed numeri- 
cally by bvp4-solver in MATLAB, by shooting in the y -P arame t er with conditions 

(2.5) f'(y ) = /'"(yo) = 0. 

The correct choice of the interface location y > is obtained from the zero-height con- 
dition /(yo) — 0, so in the limit, at y , all three free-boundary conditions (11.61) are valid. 
For numerics, we use the regularization in the fourth-order operator in (12.3ft by replacing 

(2.6) |/|" _>( e a + 

with typically e = 10~ 3 , and similar tolerances of the method. Notice a big "almost flat" 
part of the profile fo(y) in (c) for p = 2 and n = 0.8. Here, unlike other three cases, p = 2 
is sufficiently close for n + 1 = 1.8, which is another critical exponent for the ODE ( 12.3ft . 
since (5 = then. 

In Figure El we explain further details and the results of the actual shooting procedure 
for the last two previous patterns. For the reason of comparison, in (b) by dashed line we 
show the FBP profile for the semilinear case n = that was studied in [23] in the case of 
the Cauchy (not an FBP) setting. It is curious that interfaces for n = and n = 1 are 
close to each other, but not the profiles. It is worth observing that Figure El(a) reveals 
some oscillatory character of typical solutions near interfaces. This indicates that there 
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(a) p = 2, n = 0.4, y « 6.606 (b) p = 2, n = 0, y « 5.109 




(c) p = 2, n = 0.8, y « 14.822 (d) p = 3, n = 1, y « 4.455 



Figure 1. FBP profiles satisfying (12.31) . (12.41) for various values of p and n. 

exist other patterns fi for the FBP, which admit a few oscillations near interfaces and 
form in the limit the solution of the Cauchy problem; cf. [211 Prop. 5.1]. 

Let us mention again that existence and multiplicity of VSS profiles for (12. 3p are open 
problems. In the next section, we will discuss the question of existence of a finite number 
of profiles in the subcritical range in the case of the CP. 

3. VSSs in the Cauchy problem for p e (n+ l,p ) 

3.1. On local oscillatory structure near interfaces. The VSSs take the same self- 
similar form (12.11) . where the radial rescaled profile / of changing sign solves the ODE 
(12. 2p in M. N . We refer to [2"Tj § 5] and [H] for details on the oscillatory structure of such 
similarity profiles close to finite interfaces, and also to [231 Ch. 3], where "homotopy"-like 
approaches to the Cauchy problem for various TFEs are presented. 

Namely, it was shown that the asymptotic behaviour of f(y) satisfying equation (12.31) 
near the interface point as y — > > is given by the expansion 

(3-1) f(y) = (y -yY<f(s), s = ln(y - y) , fi = f , 



VSS profile for FBP 

f(y) J ' ' ' ' ^ 

0.8 - \ 

0.6 - 

0.4 - 

0.2 - 



-0.2 ' 1 1 1 1 ' 

2 4 6 8 10 

y 

(a) p = 2, n = 0.8, y = 14.822... (b) p = 3, n = 1, y = 4.455... 

FIGURE 2. Results of shooting of the profiles in Figure Q] (c) and (d). 

where, after scaling ip \— > /3nip } the oscillatory component ip satisfies the following au- 
tonomous ODE, where exponentially small as s — > — oo terms are omitted: 

(3.2) cp'" + 3(> - l)cp" + (3/i 2 - 6fj, + 2)<pt + - 1)0 - 2)y> + ^ = 0. 
According to this singularity analysis, for n G (0,nh), where 

n h = 1.759... 

is the heteroclinic bifurcation point for the ODE (13.21) . there exists a stable (as s — > +oo; 
at s — > — oo all solutions are unstable in view of shifting in yo) changing sign periodic 
solution ip(s) of (13.21) . According to (13.11) . this gives similarity profiles of changing sign, 
which being extended by f(y) = for y > y forms a compactly supported solution 
/ G C a in a neighbourhood of y = y , with a ~ -. Notice that a — > +oo as n — ► + , 
so the regularity at y = y improves to C°° at n = + . Fhese functions are oscillatory 
near the interfaces. The first results on the oscillatory behaviour of similarity profiles for 
fourth-order ODEs related to the source-type solutions of the divergent parabolic PDE 

(3.3) u t = -{\u\ m ~ x u) xxxx {m > 1), 

were obtained in [5]. It turns out that these results can be applied to the rescaled TFE 
(11.101) . but for n G (0, 1) only (the ODEs for (jl.8p and (13.31) then coincide after change). 
Some existence and multiplicity of periodic solutions of ODEs such as (13.21) are known [T5] . 
and often lead to a number of open problems. Therefore, numerical and some analytic 
evidence remain key, especially for sixth and higher-order TFEs, [23l Ch. 3]. 

Thus, we express the above results as follows: there exists a 2D bundle of asymptotic 
solutions of (12.31) in K. having the expansion (13. ip . 

(3.4) f(y) = (y — y)~tp(s + s ), with two parameters yo > and So G R, 
where we also take into account the phase shift so of the periodic orbit <p(s). 
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Therefore, matching the 2D bundle (13.41) with two symmetry conditions at the origin in 
(12. 3p leads to a reasonable well-posed problem of 2D matching, which remains essentially 
open still. In the case of analytic dependence on parameters involved, such a problem 
cannot possess more than countable set of isolated solutions. Actually, we are going to 
show that the number of VSS profiles is always finite. 

3.2. Global behaviour of VSS profiles. In Figure [3], we present the first even VSS 
profile fo(y) (and some /^(l/)) satisfying (12. 2h for N = 1, p = 2, or 3 and various n. Here 
we fix the same regularization (I2.6P where e and Tols are about 10~ 2 , which is sufficient 
accuracy. Notice that in ID, 

fo, f2, fi, ■■■ are even functions, and 
fi, h, /s, ■■■ are odd. 

As we have mentioned, the case n = corresponds to smooth VSSs for the semilinear 
parabolic equation (11. 13)) . which were studied in [24] . so we can always compare the results 
with those for the TFE (II. 3p with small n > 0. For n = 0, the VSS profiles {/;, I > 0} are 
known to appear at subcritical (for p < pi) pitchfork bifurcations at critical exponents 

(3.5) Pl = l + j^, 1 = 0,1,2,.... 

Figures |3Ja) and (b) show a strong similarity of the corresponding VSS profiles /o f° r 
various n, and, moreover, confirm that the profiles can be continuously deformed to each 
other as n — > + . This is related to a general homotopy approach to the Cauchy problem 
for TFEs and other degenerate PDEs with non-monotone operators; see [231 Ch. 3] and 

(Tsiinj. 

In Figure [3](a) for p = 2, we also show the smaller VSS profile f2{y) for n = (the 
dashed line), and also a couple of profiles for negative n = —0.1 and —0.2. These are 
/ 2 (l/), and not fo(y). In (b), we also calculate f for the negative n = —0.2. Recall that 
for n < 0, (11.31) demonstrates typical features of a fast- diffusion problem. There is no 
finite propagation in this case, but solutions are equally oscillatory as y — > +oo, thus 
inheriting this property from n = 0. Notice that for n = 0, according to (13.51) . 

p 2 = 1 + | < 3, 

so that f'2{y) does not exist for p = 3. But fi exists for positive n = 0.5 (the dotted 
line in (b)) and is sufficiently small, meaning that the corresponding critical bifurcation 
exponent p2{n) is slightly above 3. 

3.3. On a boundary layer as p — > n + 1. For n = and any p e (l,po), there exists a 
finite number 

M^^Et^+oo asp^l+ (n = 0), 

of different VSS profiles, which are obtained by standard bifurcation theory, [24] . We 
expect that a similar multiplicity property remains valid for the TFE for n > 0, though 
in this case bifurcation branches are governed by the linearised TFE operator (see [211 
§ 2] and the results below), so that bifurcation points {pi} are not given explicitly by a 
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discrete spectrum {A;} of a non self-adjoint operator as in the semilinear case n = 0; cf. 
[2U Lemma 4.1]. We analyze this kind of bifurcation in the next subsection. 

A typical strong oscillatory behaviour of the VSS profiles for p « (n + 1) + is shown 
in Figure HI for p = 2 and n = 0.95. We present here first six even VSS profiles from 
different p2/-branches; see further explanations below. Notice formation of an interesting 
"boundary layer" as p —>■ n + 1, where the VSS profiles f(y) become more and more 
oscillatory reflecting the fact that a suitable solution satisfying at p = n + 1 the ODE 
(l2~3l) for N = 1, 

(3.6) -(\f\ n ry+y-\f\ n f=Q & r, 

does not exist. Numerics in Figure 0] suggest that solutions of f !3.6j) are highly oscillatory 
and are not compactly supported. Recall that the oscillatory structure ( 13. II) near inter- 
faces demands the linear term +/3f'y in the ODE (12. 3p with (3 > 0. For (3 = 0, such 
solutions do not exist. Moreover, for /3 < 0, there exist positive solutions near interfaces, 
which correspond to blow-up problems [To] . 

3.4. Bifurcation of the first p-branch at p = p^ : a nonlinear version. We now 

study the behaviour of the first p-branch of the VSS profiles / = fo(y), when p approaches 
from below the critical exponent (11.41) . For n = 0, such a behaviour was studied in [24] 
by a standard bifurcation approach showing that the VSS profiles vanish with the rate: 

N 

(3-7) H/lloo ~ (p -p)~ -> asp^pQ. 

The bifurcation analysis in [21] was based on known point spectrum (13.191) and other 
spectral properties of a non-self-adjoint linear operator B, which is ( ll.lOp for n = 0. The 
operator B appears in the ODE (ll.lOp for n = generating the rescaled kernel of the 
fundamental solution of the bi-harmonic operator D t + A 2 . 

Local bifurcation from po- We now perform a formal nonlinear version of such a 
bifurcation (branching) analysis for n > 0. As usual, according to classic branching 
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m=2, N=1 , n=0.95, p=2: six even VSS profiles 




y 

FIGURE 4. Oscillatory behaviour of the VSS profiles for p ~ n + 1 increases: 
first six even VSS solutions for p = 2 and n = 0.95. 

theory [521157] , a justification (if any) is performed for the equivalent quasilinear integral 
equation with compact operators. For simplicity, we present basic computations for the 
differential version. 

We introduce the small parameter e = po — p, so that, as e —* 0, 

= iT^v + ^(irSvF + and = iTb - ^IPW + °(^)- 

Substituting this expansion into (12. 2p and performing the standard linearization yields 

B(/) + e(df - \f\ n ^f\n |/|) - |/r*/ + ... = 0, 
where C x = j^fiNI - f y ■ V) 

is a linear operator, and B is the rescaled operator (jl.lOp of the pure TFE. Notice that 
the fact that the operator B in (13. 8p occurs in the rescaled pure TFE correctly describes 
the essence of a "nonlinear bifurcation phenomenon" to be revealed. 
Next, we use an extra invariant scaling of operator B by setting 

(3.9) f(y) = bF{y/b n '% 

where b = b(e) > is a small parameter, (3(e) — > as e — > 0, to be determined. Substi- 
tuting (13.91) into (13. 8p and omitting higher-order terms yields 

(3.10) B(F) + eCiF - b n+ ^\F\ n+ ^ F + ... = 0. 

Finally, we perform linearization F = F + Y, where F is the "fundamental", supported 
in Bi, similarity profile of the Cauchy problem for the pure TFE, i.e., satisfying (11.101) 
for N = 1. This yields the non- homogeneous problem 

(3.11) B'(F)Y + edF - b n+ ^ \F\ n+ ^F + ... = 0. 
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Here the derivative is given by 



B'(F)Y = -V • [|F| n (f (VAF)F + VAF))] + (3y ■ VF + (3NY. 

The rest of the analysis depends on assumed good spectral properties of the linearised 
operator B'(F) . We follow the lines of a similar analysis performed for the FBP case in [2T| 
§ 2], where, in the FBP for n = 1 and F given by (II. lip , the operator B'(F) turns out to 
possess a (Friedrichs) self-adjoint extension with compact resolvent and discrete spectrum. 
Such a self-adjoint extension does not exist for the oscillatory F{y). Here we need to use 
general theory of non-self-adjoint operators; see e.g., [26] . A proper functional setting of 
this operator is more straightforward for N = 1 (and in the radial setting), where, using 
the behaviour of F(y) — > as y — > 1, it is possible to check whether the resolvent is 
compact is a suitable weighted L 2 space. In general, this is a difficult problem; see below. 

We assume that such a proper functional setting is available for B, so we deal with 
operators having solutions with "minimal" singularities at the boundary Si, where the 
operator is degenerate and singular. Namely, we find the first eigenfunction ipo with 
Ao = of B'(F). Let ipo be the corresponding first eigenfunction of the adjoint operator 
(B')*(F) defined in a natural way using the topology of the dual space L 2 (F>i) and having 
the same point spectrum (the latter is true for compact operators in a suitable space [311 
Ch. 4]). Moreover, it can be seen from the divergent form of the linearised operator B'(F) 
that, after bi-orthonormalisation, 

the corresponding first eigenfunction of (B')*(F) can be taken as 

(3.12) r (y) = i- 

This simplifies the rest of computations, though one can restore these similarly for arbi- 
trary ip*, as suggested later on for finding other critical bifurcation points {pi}- 

Further, we assume that there exists the orthogonal subspace Span{^;,/ > lj-L^o of 
eigenf unctions of B'(F), and we look for solutions of (13.111) in the form 

y = cy> + w, 

with a constant C = C(e) and a function wl-tpo, i-e., (w,iJ)q) = 0. Recall again that, 
in doing so, we need to transform (13.111) into an equivalent integral equation with com- 
pact operators, but for convenience, we continue our computations using the differential 
version; see some details in [2H § 3]. 

Thus, multiplying (13.111) by ip^ in L 2 (E>i) and integrating by parts in the differen- 
tial term y ■ VF in £\F, we obtain the following orthogonality condition of solvability 
(Lyapunov-Schmidt's branching equation [37J § 27]): 

(3-13) eiC&N)fF = V**J\F\<+*F. 

Therefore, the parameter b(e) in (13. 9p for p m p$ is given by (cf. (13.71) for n = 0) 

(3.14) 6 ( £ ) = [7o ( po _ p )]l^v ? with 7o= _g_J F /J| F |n + ^ F; 
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m=2, N=1,n=0: p-branch of f.(y), p n =5 




1 1.5 2 2.5 3 3.5 4 4.5 5 



P 

FIGURE 5. N = 1 and n = 0: the first p-bifurcation branch that is originated 
at the first critical exponent po = 5~ and blows up as p — > 1 + . 



provided that J \F\ u+ nF > (not an easy inequality that can be checked numerically). 

For n = 0, a rigorous justification of this bifurcation analysis can be found in [2U 
§ 6], where a countable number of p-branches originated at bifurcation points 03. 5p was 
detected on the basis of known spectral properties of the corresponding linear operator 
(I3.18I) ; see details in [12]. For n > 0, the justification needs spectral properties of the 
linearised operator B'(F) and the corresponding adjoint one (B'(F))*, which remain an 
open problem. In particular, it would be important to know that the bi-orthonormal 
eigenfunction subset {ipi} of the operator B'(F) is complete and closed in a weighted 
L 2 -space (for n = 0, such results are available [12] )• We expect that for n ~ 0, there 
exist critical exponents for the TFE with absorption that are close to those in (13. 5p at 
n = 0. This can be checked by standard branching-type calculus; see Appendix in |20j . 
where nonlinear eigenfunctions of the rescaled PME in ~R N were studied by a branching 
approach. 

On global extension of p-branches. For comparison, we begin with Figure that 
presents the first monotone branch of VSS profiles fo(y) in the semilinear case n = that 
exists for all p G (l,Po = 1 + jj)', cf. [24|. On the vertical axis, we put ||/||oo that, for 
such profiles, is simply /(0). 

Figure M show that such monotone branches persist until n = 0.11 (a), while, for a 
slightly larger n = 0.13, we first observe a non-monotone branch of patterns fo(y) (b), 
and this persists for most of larger n's. This is a quasilinear phenomenon to be discussed 
in greater detail below. Moreover, we observe a typical turning point of the p-branch at 
p = 2.45..., which, as usual, is characterized by existence of a non-trivial centre subspace 
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of the linearized operator, 

E a(B'(f)). 

These global bifurcation diagrams are calculated with the enhanced Tols = 10~ 4 and 
small step sizes Ap ~ 10~ 3 . We claim that the turning, saddle-node bifurcation for N = 1 
occurs above the critical exponent 

n s _ n « 0.12. 

Figure [7] illustrates the vanishing behaviour of the VSS profiles fo(y) as p — > p = n + 5 
for V = 1 and n = 1. Notice that, according to (I3.14p . the rate of decay is fast, 

ll/lloo ~ (6-p)* as p->6~. 

In Figure El we show the corresponding first p-bifurcation branch that is originated at 
p = 6~. Such a behaviour of this p-branch is similar to that for the semilinear parabolic 
equation for n = 0; cf. [24J. The global behaviour of this branch is unusual: the branch 
exhibits a definite non-monotonicity and turning for p « 3.46... . 

Such non-monotonicity branching phenomena were consistent in numerical experiments. 
Similar features are shown for n = 0.5 in Figure [9j Let us discuss possible consequences 
of such a behaviour that was not available for n — 0, and hence exists in a strongly 
quasilinear case n > sufficiently large (at least for n > 0.11 as Figure [6] suggests). It 
follows from principles of general branching theory [37] that if such a p-branches vanishes, 
it must end up at bifurcation points only. 

On the other hand, we know that the profiles fo persist until the critical value p = n+1; 
see Figure HI where fo(y) is available for n = 0.95 and p = 2>n + l = 1.95. Therefore, 
if a po-branch disappear at some bifurcation point pk > n + 1, it must appear at another 
(possibly, saddle-node) bifurcation point pi < pk- As we will explain, at standard pitchfork 
bifurcations from at p = pi with I = 1,2, other types of VSS profiles fi appear, so 
that the new appearance of fo can be associated with new bifurcations and branching. 
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Therefore, fo must appear at some subcritical saddle-node bifurcation, most probably 
embracing branches of profiles /o and fi that have a similar geometric structure. 

We now answer the question posed above. In Figure [H)] we show the po,2-branch for 
n = 1 in a neighbourhood of the bifurcation point p 2 = 3.333... and of the turning point 
shown in Figure [8] This branch has two turning points which are saddle-node bifurcations. 
It follows that the even profiles fo(y) and /^(y) belong to the same po,2-branch, i.e., can be 
continuously deformed to each other as solutions of the ODE (12. 3p . Note that the pitchfork 
bifurcation at p = p 2 is now supercritical (the branch is originated for p > p 2 = 3.333...). 
Figure [IT] shows three different "/q" profiles for p = 3.6, which is shown by the vertical 
dashed line in Figure [TU1 The smallest profile is actually f2{y), which thus is homotopic 
equivalent to fo (i.e., admits a homotopic path via a family of operators). 

We expect that such saddle-node bifurcations can occur on other ^/-branches creating 
necessary profiles in different p-intervals. For instance, we observed an evidence that the 
next pi^-branch occurs. 

Thus, the po,2-branch is a closed curve. This type of closed /i-bifurcation branches were 
earlier found in [2^1 § 6.4] for VSSs with another type of parameterization in the ODEs 
like (12.31) . n = (so that (3 — j), with the change 

i f'y ^ ^f'Vi with parameter /i G (0, 7]. 

Then a /z-branch was shown to appear at a pitchfork bifurcation point and was continued 
until another, smaller one, i.e., a global continuation of such branches up to fi = | + was 
impossible. In addition, essentially non-monotone bifurcation branches were detected 
[9J p. 1802] in rather similar fourth-order ODEs associated with blow-up in higher-order 
reaction-diffusion equations such as 

Ut = —u xxxx + |u| p_1 u in Ex (0,T). 

On other p-bifurcation branches. A similar "nonlinear" bifurcation analysis can be 
performed on the basis of any suitable similarity profiles Fi(y) of the pure TFE (11.81) . 
Namely, this profile appears in the similarity solution of (jl.8p . 

(3.15) u l (x,t) = t-°*F l (y), y = x/t*, where A = ^ > 0, 

and ai G (0, -) is a parameter. Instead of (jl.lOp . the function Fi solves the following 
elliptic equation: 

(3.16) Bi(F) = -V ■ (\F\ n VAF) + (3 t VF -y + aiF = in R N . 

According to the principles of self-similarity of the second kind (a notion introduced by 
Ya.B. Zel'dovich in 1956 [39]), the acceptable values of the parameter (3i are chosen from 
the solvability of the problem (I3.16p . i.e., from existence of a nontrivial compactly sup- 
ported solution Fi(y) in ~R N . We expect existence of a countable (up to obvious scaling) 
set of such solutions $ = {ai, Fi(y)}, with most of them not being radially symmetric. 
In a certain sense, this looks like a nonlinear extension of the linear eigenvalue problem 
that occurs for n = 0, where fii = \ and (I3.16P reads 

(3.17) BiF = -A 2 F + \VF-y + aiF = in R N , 
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m=2, N=1 , n=1 : VSSs as p tends to p n =6 




FIGURE 7. N = n = 1: VSS profiles vanish as p — > p = 6. 



m=2, N=1, n=1: p-branch of f Q (y), p =6 




4.5 5 
P 

FIGURE 8. N = n = 1: the first p-bifurcation branch of profiles /o originated 
at the first critical exponent p = 6~ ; see Figure [10] for the enlarged area of the 
turning point. 



so that \i = y ~ a i are eigenvalues of the non self-adjoint operator 



(3.18) B = -A 2 + i 2/ -V + f I. 
This spectrum is discrete and is given by [12j 

(3.19) A, = 



- 7=012 

4) L — u ; L t ^1 ■■■ 1 
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m=2, N=1, n=0.5: p-branch of f Q (y), P =5.5 
0.65 1 1 1 1 1 




P 

FIGURE 9. N = 1, n = 0.5: the first p-bifurcation branch originated at the first 
critical exponent p = 5.5~ . 
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FIGURE 10. The enlarged non-monotonicity part of the bifurcation po,2-diagram 
for the ODE (|2.3p in R for n = 1. The behaviour close to p2 = 3.333... with two 
turning (saddle-node) points. 

and each eigenvalue has finite multiplicity. This determines all possible values of the 
parameters of self-similarity 

(3.20) a, = ^, Z = 0, 1,2... . 

16 



m=2, N=1, n=1: three VSS profiles 

0.7 i , 1 , , ,— 




-0.2 ' 1 1 1 1 1 1 

1 2 3 4 5 6 7 



y 

FIGURE 11. Three /o,2(y) profiles from Figure [10] for p = 3.6. 

For small n > 0, the linear eigenvalue problem ( 13. 17ft can predict the nonlinear eigenfunc- 
tions of (I3.16P by the n-branching approach in the lines similar to that in Appendix A in 
|20j, though the justification in the fourth-order case is much more difficult. 

For larger n > 0, the nonlinear eigenfunctions {i 7 }} are unknown even in the ODE case 
N = 1 or in radial setting in M . Recall that, even in the simplest case / = 1 meaning 
the 1-dipole solution F\(y), the existence of such a profile for not small n is still unclear 
mathematically; see references and results in [7] and more recent paper [8]. 

Finally, once the nonlinear eigenfunction subset {F{\ of the pure TFE (11.81) is known 
(say, for small n > 0), for each function Fi(y) supported in Bi after rescaling, one can 
develop the above bifurcation approach taking in (13.91) F = Fi for any / = 0, 1, 2, ... . This 
gives the sequence of critical exponents 

(3-21) p = B^l = ft = l=p. =► W = l + £- 

For n — 0, (I3.2ip and (13.201) yield precisely the known critical bifurcation points (13.51) 
studied in [23] . 

Therefore, a more rigorous bifurcation analysis is available for p ~ p/ and n ~ 0, so this 
means the branching approach with the parameter // = (p, n) T ~ (1 + j^tj, 0) T in M 2 . 

3.5. On other VSS profiles. As we have mentioned, for n — 0, for any p in the subcrit- 
ical range p G (l,po)j there exists a finite number of VSS profiles {fi}, and each p-branch 
is originated at bifurcation points (13.51) ; see [2U § 6]. We claim that similar VSS patterns 
also exist for n > and each fi + i(y) has a more complicated form than the previous one 
fi(y), i.e., has "more oscillatory" structure with more non-exponentially small oscillations 
and "essential" zeros. Figure [T21 shows two VSS profiles f and fx for n = |, p — 2 and 
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m=2, N=1 , n=0.5, p=2: two VSS profiles m=2, N=1 , n=1, p=2.5: two VSS profiles 




02468 10 12 02468 10 



y y 

(a)n=§, p = 2 (b)n = l, p = 2.5 

FIGURE 12. First two VSS profiles for n = §, N = 1, p = 2 (a) and for n = 1, 
AT = 1, p = 2.5 (b). 

n = 1, p — 2.5 in the ID case (they look similarly emphasizing continuity in n). It is 
difficult to check numerically whether other profiles exist. 

For comparison, in Figure d3j we present five VSS profiles for the semilinear PDE for 
n — 0, with N — 1 and p = 1.7, which are computed much easier. According to (|3.5p . 
there exist five bifurcation points above 1.7, 

(3.22) p = 5, pr = 3, p 2 = f , p 3 = 2, p 4 = | = 1.8 (p 5 = § = 1.66... < 1.7). 

Therefore, for p = 1.7, bearing in mind that the p-bifurcation branches are monotone 
decreasing in p, there exist precisely five VSS profiles indicated in this Figure. By conti- 
nuity, we expect that these profiles do not essentially change and exist for all small enough 
n > 0, where these become compactly supported with the interface dependence [T5] § 10] 

yo(n) ~ rr'~± — > +oo as n — >• + . 
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